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Motivation 


Second-order methods are currently popular; for turbulent and 
unsteady flows, they are generally unreliable or impractical. 


Worldwide effort to improve accuracy and efficiency for such 
flow problems by employing high-order methods: 


|) The numerous papers at this conference. 
2) International Workshop on High-Order CFD Methods 


3) The TILDA project (Towards Industrial LES and DNS for 
Aeronautics) supported by the European Union. 


High-order methods need further development and 
improvement. The current work is along this direction. 


Setup 


Goal: solve the Navier-Stokes equations numerically. 
Numerical methods toward this goal are typically first 
derived and analyzed for the advection equation. 


Van Leer (1977) introduced five schemes for advection in 
‘Towards the ultimate conservative difference scheme, IV”’. 


Scheme I: least accurate but became most popular; widely 
known as the MUSCL scheme. 


Schemes III (piecewise linear, discontinuous) and V 
(piecewise parabolic, continuous): most accurate but least 
popular. They are the main focus of current work. 


Setup 


(Van Leer and Nomura 2005): “When trying to extend 
these schemes beyond advection, viz., to a nonlinear 
hyperbolic system like the Euler equations, the first author 
ran into insuperable difficulties because the exact shift 
operator no longer applies, and he abandoned the idea”. 


Scheme III was extended to systems of equations in 
(Huynh 2006). The approach was further analyzed and 
applied to hyperbolic-relaxation equations by Suzuk1, 
Khieu, Van Leer (2007-2009). High-order extension was 
carried out by Lo (2011) and Huynh (2013). 


Scheme V is being extended to systems of equations by 
Roe, Eymann, and Fan (2013-present) and called the 
active flux scheme. 


Main Findings of Current Work 


Equivalence result: Schemes III and V are shown to be 
equivalent in the sense that they yield identical solutions. 


© This equivalence is counter intuitive. 


© This finding also shows a key connection between the 
approaches of discontinuous (scheme III) and 
continuous (scheme V) polynomial approximations. 


High-order extension: introduce a projection-interpolation 
framework that simultaneously extends schemes HI and V. 


Outline 


Review MUSCL approach and schemes III and V. 


A new result: the equivalence of these two 
schemes. 


Introduce a projection-interpolation framework 
that simultaneously extends schemes III and V to 
arbitrary-order. 


Von Neumann (or Fourier) stability and accuracy 
analysis. 


Conclusions and discussion. 


Advection Equation 
u.tau,=0, a2zO0 
Initial condition: at t=O, u(x) =Uuo(x) 


Exact solution: u(x,t) =up(x—at) 
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Advection Equation 


Simplifies the derivation and description of CFD 
methods. 


Facilitates linear stability and accuracy analysis. 


Schemes derived for advection must then be 
extended to systems of equations, which is often 
not an easy task. 


Discretization 


The data u j represents the average value of u in the cell 
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First-Order Upwind Method for Advection: 
Shift Operator and Projection 
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High-Order Extension: Legendre Polynomials 


On / =[-1,1], by orthgonalizing the monomials ¢€ ' ,k =0,1, 2, ..., 


via Gram - Schmidt process, we obtain the Legendre polynomials 
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Projection using Legendre Polynomials 


On / =[—1, 1], approximate u by aN jk Le where 
1,0 =f Qué |g? = 4f) tuae. 
Wi =f, bude /a/. 
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Van Leer’s Scheme III Employing Shift and Projection 
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(b) Linear solution 


Van Leer’s Scheme III 
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Scheme III 
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Set 


Then 


Projection to Arbitrary Order 
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(a) Data (b) Solution (red curve) 


Scheme III can be considered as a precursor of the DG method. 


Van Leer’s Scheme V 


The cell average values u 7.0 and interface values u j+1/2 are stored. 


In each cell 7, the values u j-1/2? 4;,0° and u j+1/2 define a parabola. 
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(a) Parabolic Data (b) Parabolic solution (red) 
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Equivalence of Schemes III and V 


* Assume that the cell average values u 1,0 and interface values u j+1/2 
for scheme V are known. 
* Also assume that the CFL number o 1s fixed. 
* For scheme III, set 
Uj, = A-O)U 5,172 “Uj 9) + TU; 9 —45_1/2) 
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* Then the cell average solutions at time t’" of schemes III and V are 


identical : 
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Equivalence of Schemes III and V 
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Equivalence of Schemes III and V 


* Assume that the cell average values u 1,0 and interface values u j+1/2 
for scheme V are known. 
* Also assume that the CFL number o 1s fixed. 
* For scheme III, set 
Us = A-O)U 54472 Uj 9) + CU; 9 —4j-1/): 


n+l 


* Then the cell average solutions at time t’" of schemes III and V are 


identical : 


a 
uid — uig= Uso: 
* In addition, 
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Projection-Interpolation Schemes 


| piece of data per cell,or K =1 
n+1 
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Projection-Interpolation Schemes 


2 pieces of data per cell,or K =2 


Projection P, Poly I, or Cubic 
Scheme III Scheme V Interpolation 
» e 


Projection-Interpolation Schemes 


3 pieces of data per cell,or K =3 
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Family of PuIv Schemes 


* Projection part: Projecting the data onto P,, resultsin 44 +1 quantities 
Wi =Us x, R=0, ni 


where 


2k +1 
a — 5 (u, L, ). 


* Interpolation part: At each interface x; ,;;2, the interpolation part 
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consists of approximations to = degree v,1.e., the v+1 quantities 
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Family of Pulv Schemes 


* For each cell, we have w+1 projection quantities 
in the celland 2(v+1) interpolation quantities at 
the two interfaces. 

* These 2v+ +3 quantities for each cell define a 
polynomial of degree 2v+ w+ 2. 

* Shift the polynomial data a distance aAt. 

* Update the projection part and the interpolation part. 
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Fourier Stability and Accuracy Analysis 


Define U; with K =4+v+2 components by 
U. = 


J (Uso Wh rra0™ serra vl 
The solution can be written as 
n+l 

Assume the data is a harmonic that satisfies 

lw 
U iy = 2 U 
Then, 

um! - eC, + e!C 14+ CU; 
For the first - order upwind scheme, 


Une = Pau 4 d—o)|U,. 
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Fourier Analysis: Plots of Absolute Values of Eigenvalues 


Fourier Analysis: Plots of Absolute Values of Eigenvalues 


Fourier Analysis: Plots of Absolute Values of Eigenvalues 
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Projection-Interpolation Schemes 


Derived for advection equation 


With K degrees of freedom per cell, the scheme is 
accurate to order 2K — 1 (..e., it is Super accurate or 
Super convergent) 


CFL condition is 1 as opposed to ~1/K* for explicit 
RK-DG. 


Extension to systems of equations in multiple 
dimensions remains an open problem. 
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Conclusions 


Reviewed MUSCL approach for the advection equation. 


Presented a key new result: the equivalence of schemes III 
and V. 


The above result shows a key connection between 
continuous and discontinuous approaches. 


Introduced a projection-interpolation framework that 
simultaneously extends these schemes into a single family 
of high-order methods. 


Discussed Von Neumann (or Fourier) stability and 
accuracy analysis. 


Further research on these methods is needed. 


Thank you 
for your attention. 


Examples for the Equivalence, o =0.2 
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Examples for the Equivalence, o =0.8 
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